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Modern applications require a robust and theoretically strong tool for the realistic modeling of 
electronic states in low dimensional nanostructures. The k ■ p theory has fruitfully served this role 
for the long time since its creation. During last two decades several problems have been detected in 
connection with the application of the k ■ p approach to such nanostructures. These problems are 
closely related to the violation of the ellipticity conditions for the underlying system, the fact that 
until recently has been largely overlooked. We demonstrate that in many cases the models derived 
by a formal application of the Luttinger-Kohn theory fail to satisfy the ellipticity requirements. The 
detailed analysis, presented here on an example of the 6x6 Hamiltonians, shows that this failure has 
a strong impact on the physically important properties conventionally studied with these models. 
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The effective mass theory is one of the fundamen- 
tal parts in the physics of nanostructures. This the- 
ory allows us to get theoretical insight into the elec- 
tronic properties of dominating bands near the extremum 
point^l^. Futhermore, the theory establishes a ro- 
bust computational framework for simulating observable 
quantum-mechanical states and corresponding energies 
in the low-dimensional systems including quantum wells, 
wires, nanodots. In the original Luttinger-Kohn worls^l 
authors applied the theory to the Schrodinger equation 
perturbed by a smooth potential and constructed a repre- 
sentation for valence bands Hamiltonian near high sym- 
metry point F of the first Brillouin zone in bulk Zinc- 
Blende (ZB) crystals with large fundamental band gap. 
Soon after that, Kane showed how to expand the model 
to the narrow gap materials such as InSb and Ge for in- 
stance, where one can also account for the influence of 
the conduction band^. One of the advantages of the 
k ■ p theory is in its universality. Indeed, the theory 
had also been extended to cover Wurtzite (WZ) type of 
crystals, materials with inclusions, heterostructure ma- 
terials and superlatice^. Another advantage of the ef- 
fective mass theory that has recently been explored in 
some details in IHl is its flexibility, as one can easily ad- 
just the models to include additional effects like strain, 
piezoelectricity, magnetic field, and respective nonlinear 
effects. These inbuilt multi-scale effects are crucial for 
such applications as light-emission diodes, lasers, high 
precision sensors, photo-galvanic elements, hybrid bio- 
nanodevices, and many other^. 

For a wide range of applications these models 
have provided good, computationally feasible and effi- 
cient approximations that agree well with experimen- 
tal results^ ^. However, for some types of crystal ma- 
terials band structure calculations based on such multi- 
band models lead to the solutions with unphysical 
properties^ or so called spurious solutiond^^J^i^J 

As a result, there have been various attempts to ex- 
plain the origin of the spurious solutions and develop 
some reliable procedures on how to avoid thenP. These 
approaches rely on three main ideas: (a) to modify the 
original Hamiltonian and remove the terms responsi- 



ble for the spurious solut ion j^^*^, (b) to change band- 
structure parameter j^^ l ^^ l, and (c) to identify and exclude 
physically inadequate observable stateJ^^. All mentioned 
approaches suffer from the common weakness - the lack 
of clear justification of the underlying theoretical proce- 
dure and thus from limitations in their applicabilitjISIIll 
In this work we show that spurious solutions are not 
just a reason but rather an imminent consequence of an- 
other fundamental problem in applications of the classi- 
cal theory - the non-ellipticity of the multiband Hamilto- 
nian. The widely adopted effective mass approximations 
of the original Schrodinger elliptic Hamiltonian turns out 
to be non-elliptic for a broad class of known material pa- 
rameters (cf. Table |l]). This fact leads to the following 
consequence: since any qualitative approximation meth- 
ods must preserve the topological structure of the spec- 
trum of a general linear operator, and as an implication 
symmetric properties in case of Schrodinger Hamiltonian, 
this evident discrepancy signifies the mathematical inva- 
lidity of performed approximation procedures for those 
materials. Among these mutiband Hamiltonians there 
are a number of widely used 6x6 and 8x8 Hamiltoni- 
ans. 

The paper is organized as follows. First, we revise 
basic properties of the original Schrodinger equation, and 
outline the mathematical model of the electronic band- 
structure problem. 

Next, we derive the exact constraints on the material 
parameters for typical 6x6 Hamiltonians in ZB"^ * and 
materials. Direct calculations using conventional 
Luttinger parameters (e.g. I19I20I) show that such pa- 
rameters (e.g. Table |l]) for many important semiconduc- 
tor materials entail the violation of ellipticity require- 
ments. As a result, the corresponding bandstructure 
model for them potentially susceptible to unphysical so- 
lutions, even in the bulk case, and therefore ought to be 
modified. Moreover, the corresponding time dependent 
Schrodinger equation loses the fundamental property of 
state conservatiorPiJ. 

The material properties (such as fundamental band- 
gaps and spin-orbit splitting energies) obtained exper- 
imentally, represent real phenomena, whereas models 
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based on finite bands Hamiltonians are meant to approx- 
imate them. The very last step in such approximation 
scheme^i^ enables us to calculate interband corrections 
to the main par t of th e Hamiltonians with help of per- 
turbation theorjllJ2122] However, this last step lacks a 
rigorous theoretical foundation as it does not guaranty 
the convergence of the perturbation expansion. The re- 
sult is that the derived Hamiltonian, although directly 
based on experimental parameters (column 4, |l|, repre- 
sents a totally different mathematical object compared to 
its origin. The physical evidences, to support this claim 
have been already known for GaA^2ll and recently been 
reported for SP^. 

We start with the Schrodinger equation with the po- 
tential V{x) 

Hoi>ix) = ^i^i^) + iV{x) + Eimx) = Ei^ix), (1) 

where is the band edge energy of the system, p 
is a momentum operator, x € C C",n < 3, m* 
is a piecewise-constant effective mass parameter. In f2 
we supplement ([T]) by usual Dirichlet conditions on the 
boundary S = dft 

^{x) = 0, x(zS. (2) 

Our major focus is on the heterostructure case, where 
the difference in the effective mass and probability cur- 
rent conservation imposes the discontinuity of V?/'(x) on 
the interface between materials (assuming that the basis 
Bloch functions are the same in all constituent J^Il) jn 
this case an appropriate choice of the functional space 
for coefficients from ([I]) is the Hilbert space of distribu- 
tions with the compact support •7^i(J7pfil22I^ 

Then, the operator Hq with the domain of definition 
D{Hq) C ^{^{fl) is a symmetric operator with an existing 
self-adjoint extension, hence it conserves the probability 
currenli2Si. General theory of elliptic partial differential 



equations characterizes the problem ([T]) in terms of the 
corresponding Fredholm type theorems (e.g. [29!), so that 
the problem Q, ([2| has a countable set of eigenvalues 
Ei, smallest eigenvalue -Bo is simple and corresponding 
eigenstate "00 is of constant sign in fi, furthermore 

Ei£[pQ,oo), pq = E^ + Tmn{0,miV{x)}, i = 0,...n, 

where V(x) is real. 

If V{x) is a gently varying function over the unit cell 
in the sense of |3| the original operator Hq can be approx- 
imated by another operator H (using Bloch theorem), 
determined by the projection P of Hq on t he c onsidered 
eigenspace and Lowding perturbation theorjfflll. The last 
step in this approximation procedure accounts for the in- 
fluence of the elements from the space complement to the 
egeinspace by the formula 

r 

H = PHa + ^ yn^"^^ (3) 
1=1 

up to the order r. Setting A = 1 leads one to the final 
approximation, under the assumption that the series ([3]) 
is convergent for such A. Despite a wide applicability 
of such approximations, the intrinsic ellipticity require- 
ments for the realizations of H have not been explicitly 
verified in a systematic manner (see [T k4l221 as well as 
more recent works). 

As an example, let us consider two classical Hamil- 
tonians for ZB^ and type of materials with more 
scrutiny. In what follows we use the parameter notation 
identical to the works where corresponding Hamiltoni- 
ans were obtained. When necessary, the parameters will 
be converted from one notation to another by using the 
formulas from [p. 82. [50]. 

First consider the Luttinger-Kohn (LK) Hamiltonian 
from [3] 
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where each of the P,Q,L,M is a second order position 
dependent differential operator or equivalently second or- 



I 

der polynomial in the momentum representatiorP: 

P = {A + B)k+k^ +2Bkl, Q = Bk+k^+Akl, 

= ^((^ - ^)(^' " - 2zCfc,^,), 
L — i 2^/3 ^— , ^± — kx ^ iky. 
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Our aim is to check the type (elhptic, hyperbohc or es- 
sentially hyperbolic) of the H^^ as a partial-differential 
operator (PDO) on 'H^{{1) as we know that the given 
Schrodinger operator from ([T]) is elliptic. Only the sec- 
ond order derivative terms are playing the dominant 
role in the following analysis because contributions from 
the terms linear in k as well as from the potential, are 
bounded in the domain D{H^^j^. It means that the 
results for more complicated physical models with po- 
tential contributions from additional fields (e.g. strain, 
magnetic field, etc.) will stay the same as for the original 
ff^^, analyzed here. The fact that the Hamiltonian is 
a linear operator guarantees that it is also true for any 
other representation of H^^ obtained by linear (basis) 
transformations. 

In a more general sense, for any m-dimensional matrix 
PDO H — {hijYl^j^i, where each element hij is a second 
order one dimensional PDCJ^SEH 

k,l=0 

the associated quadratic form is defined by 

GiCu-,Cnm) =vMv'^, V=i^l,...,^nm), (5) 

where M is an mn x mn matrix composed from the el- 
ements h'fj . The k • p Hamiltonians in M'^ are a special 
case of (|4])-([5]) with the real eigenvalues and n = 3 (e. 

g-m- 

Using these notations, the procedure of obtaining the 
ellipticity condition for H reduces to the question about 
the sign of for the associated M. More precisely, the 
matrix differential operator H will be elliptic if and only 
if all eigenvalues o f the corresponding Hermitian M will 
have the same sig 

In general, it is a challenging task to calculate the 
eigenvalues of M explicitly, even for such small as 3 x 3 
Hamiltonians, but we recall that we deal here with usu- 
ally sparse, band structure operators. 

Taking into account the fact that the sequence of 
eigenenergies of Hq is semi-bounded, for an approxima- 
tion H^^ (in the momentum representation), we obtain 

< 0, Vi ==0,l,...,nm. (6) 

The last constrains guarantee the ellipticity (in strong 
sens^^H) of Hamiltonian H. The operator H possess a 
self-adjoint extension in D{H) C H''^'^{fl), fc > 0, pro- 
vided that the domain is sufficiently smooth (piece- 
wise Lipschitz). Then it can be extended to a Hermi- 
tian operator by closure in the norm [p. 113, 29 or via 
Lax-Miligram procedur^^ll From the physical point of 
view the smoothness characteristics of D{H) fulfill the 
natural assumption of quantum theory that the state 
of the system must be a continuous function of spatial 
variables even when some coefficients of H have finite 
jump^^H like in t he heterostructures consisting of differ- 
ent materialJ^Ml. 



TABLE I: The material parameters for ZB type of materials, 
d - distance from the point {A, B, C) to the ellipticity region 



El. 
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B 
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d 


El. 


A 


B 


C 


d 


AlAs" 


-7.5 


-2 


8.4 


1.97 


Ge-^" 


-30 


-4.6 


33 


10.64 


AlP" 


-3.8 


-3.4 


7.2 


0.92 


Ge''" 


-30 


-4.6 


33 


10.64 


AlP^ 


-3.7 


-3.4 


7.2 


1.01 


Ge" 


-30 


-4.4 


36 


11.99 


GaN" 


-7.5 


-3.8 


6.1 


In 


InP^ 


-11 


-1.8 


10 


2.86 


CO 


-4 


-3.4 


6.6 


0.42 


InP^" 


-15 


-2.1 


17 


5.72 


C'' 


-2.9 


-2.3 


3.9 


In 


InP^ 


-9 


-3.3 


9.6 


1.34 


C'' 


-6 


-3.8 


6 


In 


InSb" 


-100 


-4 


96 


39.35 




-3.1 


-1.7 


0.9 


In 


InSb" 


-100 


-3 


96 


40.25 


GaAs' 


-16 


-2 


6 


0.89 


Si* 


-6 


-3.5 


9.6 


1.18 


GaAs' 


-17 


-2.2 


6.6 


0.98 


Si* 


-5.5 


-3.6 


8.4 


0.54 


GaAs' 


-15 


-3.1 


17 


4.83 


Si* 


-6 


-3.4 


8.4 


0.72 


GaP' 


-6 


-3 


7.2 


0.54 


SiC" 


-2.6 


-1.7 


4.2 


0.38 


GaP'" 


-8 


-2.2 


10 


2.50 


SiC^ 


-4.8 


-1.8 


5.1 


0.67 



''Measured under T = 4:.2K 

"^Obtained by extrapolations from 14 X 14 k ■ p model 

•^Measured under T = 300K 

"Ref.^O 

''Refiia 

'^RefPni 

^Most probable value (set 5 fronflSl) 
''Sets 1, 2, 3 from^i 

'Sets 7 (T=50K), 8 (T=70K), 2 fronPl 
^Sets 7 (T=60..300K), 2, 1 fronPl 
^Sets 2 (T=1.26K), 3, 6 froniSl 
'Sets 4 (T=1.6K), 3, 6 fronJSl 

The direct calculation by ([s]) for H^^ {n — 3, m — 6) 
leads us to the 18 x 18 matrix A/^^ with the following 
distinct eigenvalues: 

Xi = -§ + A, A2/3 = (|±^)C + ^, 

(7) 

•^4/5=±T + -8i •^6/7 = ±T + '^' 

where Ai, A2/3, A4/5 have the multiplicity 2 and Ag/7 
have the multiplicity 4, respectively, and A,B,C are 
usual material parameter^. By substituting ^ into 
([g]), we receive the system of linear inequalities with 
respect to A, B and C. They describe the feasibility 
region in the A, B, C space, when H^^ is an elliptic 
partial differential operator with the discrete and de- 
creasing sequence of eigenvalues. One can use similar 
reasoning to obtain corresponding inequalities for other 
common representations of through Luttinger pa- 

rameters jp^. Evidently, any solution of (|6| for Q 
would have a unique corresponding solution in 7^ no- 
tation. Note next that the solution domain A-{H^^) 
of ([7|, Q is symmetric with respect to the sign of C. 
It follows from the form of the H^^ which is depen- 
dent on |C| only. A_(7?^^) comprises an unbounded 
pyramid in M.^ (cf. Fig. [l]) with the following rays as 

its edges: {-^t,-^t,-t) , (- + ^) t, -i^, i) , and 
i-t, 0, 0), (0, - (f + ^) t, 0) , forte [0, 00). 



4 




FIG. 1: The region A_(i/^^) and the material parameters 
from Table |l] (color online) 

Figure [l] represents the H^^ ellipticity region together 
with the widely adopted values of the material parame- 
ters for different ZB type materials summarized in Table 

in 

One can observe that, among all analyzed materials 
only two (indicated as "In" in the third column of Table 
|l]) have the admissible sets of parameters. All other data 
from Table [l] yield the condition A4 > 0. That is why 
the H^^ , for the corresponding materials, is not elliptic 
and may not be even symmetric. Moreover, instead the 
domain D{H^^) = D{H) C n^ifl) we have only 

D{H^^)^ D{H) (zH\n), (8) 

it means that the solution of ([ij will have the discontinu- 
ities, for the systems with jump discontinues coefficients, 
which is the case for heterostructure materials, see e.g. 
[321 The general theory guaranties that, in this case, the 
interface discontinuities will be observed all through the 
interior of the active region along the characteristics of 
i/^^, which are now shown to exist since the associated 
G in ([5| is of nonconstant sign [p. 153, [29]. Additionally, 
([S]) and the double degeneracy of A4 ([T]) would lead to the 
nonanalytic solution and thus the momentum operator, 
from ([I]), will be ill-defined (by the embedding theorems, 
[p. 119, l29J). All the arguments stated before allow us 
to conclude that the H^^ does not provide a sufficiently 
good approximation, preserving the type of the PDO, for 
most of the practical data. 

The ellipticity analysis for ZB, can be applied to a WZ 
6x6 Hamiltonian^ without any changes. Ellipticity 
conditions that follows from such analysis are, again, lin- 



ear in parameter variables 

A2 + Ai< 0, A2 ± {2A5 - 2A4) < 0, 

(9) 

A2 + AAz < 0, A2 + -iAi - 2A5 < 0. 

Namely, Ai are well-known Luttinger-like parameters 
for WZS. As in the ZB case, each separate inequality 
has been obtained from ([6]) by substituting every dis- 
tinct eigenvalue of the matrix M associated with the 
WZ Hamiltonian. These conditions are also violated in 
most of practically important materials, among which, 
we would like to mention GaN, AIN and ZnO. Here, the 
distances to the WZ ellipticity region defined by ^ are 
approximately equ al: 0.8 04, 0.862, 0.606 for GaN param- 
eter setPHlE i.is^Ml]^ 1.27ll^, l.OllSHfor AIN; IMI^ 
for ZnO. The distances have the order of terms in the un- 
perturbed part of the WZ Hamiltonian (which in the di- 
mensionless Luttinger-like notation equal to 1), and thus 
are considerably high. 

Let us return back to the feasible parameters for two 
materials C and GaN. For carbon, the parameter val- 
ues were analyzed in SS, where authors showed that they 
don't agree with the Hall effect experimental measure- 
ments. In the same paper the authors suggested another, 
more consistent (in term of the measurements), set of pa- 
rameters (row 5, Table[l|. Observe, however, that the lat- 
ter set does not belong to the ellipticity region A_(//^^). 
In terms of the distance to A_ (H^^), we can also classify 
other mostly large band gap materials, such as Si, SiC, 
AlP and GaN, as those belonging to the same group. For 
GaN we have the set of A, B, C lying inside the region 
and for other three materials sets lie relatively close to 
this region. Such small deviations are within the reported 
order of measurement accuracy (0.25%, 0.3% 0.6% for Si, 
SiC and AlP, respectively). They can be eliminated by 
direct adjustments. The fact that the Si belongs to that 
category in spite of its smaller band gap of ~ 1.11 can be 
easily explained. Indeed, it is one component diamond 
crystal with highly regular parabolic main valence and 
conduction bands diagrams, and additionally its struc- 
ture follows the time reversal symmetry at F point. 

The rest of the materials from Table U have more com- 
plicated structure, e.g. the InSb is a small band gap, 
big effective mass material. It is knowrP, that by ac- 
counting for the valence bands only, LK approximation 
would be insufficient for InSb like materials, and pre- 
sented analysis support this fact theoretically. Concern- 
ing the Ge and GaAs they have anisotropic lower con- 
duction bands without time reversal symmetry and high 
coupling between the p-bonding topmost valence band 
and p-antibonding conduction band state^^. Inclusion 
these conduction band states leads to more precise 8x8 
and 14 x 14 models'* The setup described above is 
still applicable for these models but with a few minor 
modifications. 
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